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We study the quantum phases of anisotropic XY spin chain in presence and absence of adiabatic 
quench. A connection between geometric phase and criticality is established from the dynamical 
behaviour of the geometric phase for a quench induced quantum phase transition in a quantum spin 
chain. We predict XX criticality associated with a sequence of non-contractible geometric phases. 
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Introduction: In recent times, the quantum phase tran- 
sition (QPT) has become one of the prime research top- 
ics in condensed matter physics both from the theoret- 
ical and experimental perspective. QPT is associated 
with the fundamental changes that occur in the macro- 
scopic nature of the matter at zero temperature due to 
the variation of an external parameter. Quantum phase 
transitions are characterized by the drastic change in the 
ground state properties of the system driven by the quan- 
tum fluctuations. 

In this letter, we study the quantum phases and quantum 
phase transition of exactly solvable XY spin chain model. 
Although the exactly solvable XY model has been well 
studied [l| and is known to present a very rich structure 
but still the effect of linear quenching on the quantum 
phases has not been studied explicitly. The study of 
quantum phase transition and the nature of criticality 
through the analysis of the dynamics of geometric phase 
for linear quenching process is also rare for this model. 
The XY model exhibits different regions of criticality, like 
XX criticality, XY criticality and Ising criticality depend- 
ing on the values of the parameters (anisotropic exchange 
interaction and magnetic field) of the system. In this let- 
ter, we study the nature of criticality explicitly through 
the dynamics of geometric phase when the system under 
consideration is under quench induced QPT. We also ad- 
dress the issue on the nature of the geometric phase in 
the context of XX criticality; it is shown that the XX 
region of criticality is characterized by the existence of 
non-contractible geometric phase. Here we mention very 
briefly the essence of geometric phase in condensed mat- 
ter: Geometric phases have been associated with a vari- 
ety of condensed matter phenomena 0-0] since its incep- 
tion [§]. Besides, various theoretical investigations, geo- 
metric phases have been experimentally tested in various 
cases, e.g. with photons |9l4ll|. with neutrons 12|, |l3| 
and with atoms 1J| . The generation of a geometric phase 
(GP) is a witness of a singular point in the energy spec- 
trum that arises in all non-trivial geometric evolutions. 
In this respect, the connection of geometric phase with 



quantum phase transition (QPT) has been explored very 
recently 15-l3]- The geometric phase can be used as a 
tool to probe QPT in many body systems. Since response 
times typically diverge in the vicinity of the critical point, 
sweeping through the phase transition with a finite veloc- 
ity leads to a breakdown of adiabatic condition and gen- 
erate interesting dynamical (non-equilibrium ) effects. In 
the case of thermal phase transitions, the Kibble-Zurek 
(KZ) mechanism [11, explains the formation of de- 
fects via rapid cooling. This idea of defect formation in 
second order phase transition has been extended to zero 
temperature quantum phase transition (QPT) [13, HH by 
studying the spin models under linear quench. We will 
use this concept in our study. 

Quantum Phase Analysis and Effect of Linear Quench: 
Let us start with the model Hamiltonian 



M 



H = 



E (- 



-ai CTi+i 



(1) 

where i is the site index, x, y, and z denote components 
of spin, a is the anisotropic coupling strength and B 
is the linear quench induced magnetic field in the z 
direction. 

We recast the spinless fermions operators 
in terms of field operators by the relation 
'ip{x) = [e^kpx ^j^l^x) + e^^'^"^ i'Lix)] where 
iJjr{x) and tphix) describe the second-quantized fields of 
right- and left-moving fermions respectively. One can 
express the fermionic fields in terms of bosonic fields by 



the relation ibrix) = -H= 



(r0(x) - e(x))^ ^ denoting 

the chirality of the fermionic fields: (-fl) for right or 
(-1) for left movers. The operators Ur commute with 
the bosonic fields. Ur of different species commute and 
Ur of the same species anticommute. (p field corresponds 
to the quantum fluctuations (bosonic) of spin and 9 is 
the dual field of 0. They are related by the relation 
(j)R = 9 — (f) and 4>l = 9 + (j). We finally get the 
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bosonized Hamiltonian as 



H = Ho + —r^ / dx cos[2J — e{x)] 

+B{t)^ J dx d,^{x). (2) 



In this derivation, we have used the foUowing expres- 
sions for spin operators in terms of the bosonic fields: 
= [ C2COs(2V^0) + (-1)"C3 ] cos(yf^) , 

sin(^0), = 




5^ = C2Cos(2V7ri^0) 

d,j;4' + (-l)"ciCOs(2V7rZ0) . if is the Luttinger 
liquid (LL) parameter, C2 and C3 are the constants. 
Now we analyse the quantum phases of the system in 
the absence of linear quench: The second term of the 
Hamiltonian, which is the sine-Gordon coupling term is 
relevant when K > 1/2. The elementary excitation is 

gapped (one can estimate the mass gap by this relation 

1 

M = A(^) 2-i/-ff , A is a cutoff parameter) and the sys- 
tem is in the staggered order phase. On the otherhand, 
for K < 1/2, the system is in the Luttinger liquid phase. 
Next we interpret the quantum phases of the system by 
analyzing the two renormalization group equations. The 
renormalization group (RG) equation for the Hamilto- 
nian (Eq. 2) is 



da 



l/K)a, 



dK 
If 



aV4. 



(3) 



We interpret from the first RG equation that sine-Gordon 
coupling with strength a will be relevant and the system 
will fiow from the weak coupling phase to the strong cou- 
pling phase, when K > 1/2 otherwise this coupling term 
is irrelevant and the system is in the LL phase of the 
system. The second RG equation reveals that as a will 
increase, the LL parameter will increase, i.e, the fiow of 
the second RG equation also support the flow of the first 
RG equation for its relevant phase, i.e., the system is in 
the staggered order phase. Results derived from Abelian 
bosonization method is consistent with the RG study be- 
cause these RG equations have only trivial fixed point. 
Now we analyze the effect of linear quenching (last term 
of Eq. 1) , it modifies the quantum phases of the system. 
In the static limit for K < 1/2, the system is in the LL 
phase. The system drives to the ferromagnetic phase due 
to the presence of linear quench induced magnetic field 
in the z direction. 

In the static limit, sine-Gordon coupling term is relevant 
{K > 1/2), the system has a excitation gap (stagger 
phase). The system drives to the ferromagnetic phase 
when the quench induced magnetic field is larger than 
the excitation gap of the system. The system is in the 
LL phase when the quench induce magnetic field is of the 
order of magnitude of the gap {B{t) M = A(f )5^^). 
At time t = 0, system is in the staggered order magnetic 
phase. 



Geometric Phase and Criticality: In this model, the geo- 
metric phase of the ground state is evaluated by applying 
a rotation of 4> around the Z-eods in a closed circuit to 
each spin 1^ 23] ■ A new set of Hamiltonians is con- 
structed from the Hamiltonian (1) as 



(4) 



where U (0) — n7=-j\/ exp(i(/)cr| /2) and (t| is the z 
component of the standard Pauli matrix at site j. 
The family of Hamiltonians generated by varying 
has the same energy spectrum as the initial Hamilto- 
nian and 77(0) is 7r-periodic in 4>. With the help of 
standard Jordan- Wigner transformations, which makes 
the spins to one dimensional spinless fermions via 

the relation a.j = (ni<j '^t^ ^^'^ then using the 
Fourier transforms of the fermionic operator, dk = 
Ej exp (^f^) with k = -M, ... + M the 
Hamiltonian can be diagonalized by transforming the 
fermionic operators in momentum space and then using 
Bogoliubov transformation. The ground state \g > of the 
system is expressed as 

|5>= 11(^0^ yI'^>''^ |0>-fe -^e'^^sin^|l >, |1 >^k) 

k>0 

where |0 >k and |1 >k are the vacuum and single 
fermionic excitation of the fc-th momentum mode respec- 
tively. The angle 9k is given by 



cos fJk 



cos k — B 
A^ 



(6) 



and Kk — \J (cos k — BY + oi^ siv? k is the energy gap 
above the ground state. The ground state is a direct 
product of N spins, each lying in the two-dimensional 
Hilbert space spanned by |0 >fc |0 >_fc and |1 >k |1 >~k- 
For each value of fc, the state in each of the two di- 
mensional Hilbert space can be represented as a Bloch 
vector with coordinates {2(f>,9k). The overall phase is 
given by the sum of the individual phases. One can 
also write the above Hamiltonian as single particle ex- 



citations H{a,B{t),(j)) = ^kbk^bk, where bk = 

cos{9k/2)dk — ie^'^'^ sin{9k/2)dk'^ . The pseudomomenta k 
take half integer values: k 



27r 



AT-l 27r 



The 



"2 AT ' ' -^2 N ■ 

direct calculation shows that the geometric phase for the 
kth mode, which represents the area in the parameter 
space enclosed by the loop determined by {2(f), 9k) is given 

by 



Ffe =7r(l-cos0fe) 



(7) 



The geometric phase of the state \g > is given by Tg — 
J2k Tfc. For an adiabatic evolution, if the initial state 
is an eigenstate, the evolved state remains in the eigen- 
state. So we may now derive the instantaneous geometric 
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phases of this system due to a gradually decreasing mag- 
netic field. Let us explore the situation when the system 
(1) is driven adiabatically (slow transition ) by a time 
dependent magnetic field B{t) such that 



B{t < 0) 



t 



(8) 



B{t), driving the transition, is assumed to be linear with 
an adjustable time parameter Tq (l/r^ is the quenching 
rate). Let the system be initially at time t{< 0) << Tq 
such that B{t) >> 1. The instantaneous ground state at 
any instant t is given by 



^|0>, |0>_,-ze(2^«sin^|l>, 

(9) 

We now use eqn. (7) and (6) to derive the geometric 
phase of the fc*'' mode which yields 



cos fc + — 



(cos k + — )2 + q;2 sin k 



(10) 



The geometric phase for an isotropic system with a = 
and quantum Ising model with a = 1 may now be easily 
obtained. 



For a = 0, rk(t) = 27re(|t| - and 

cos k - 

for a = 1, rk(t) = TT I 1 - 



1 



2— cos/c 



11) 



For a system of size N the total geometric phase for the 
initial state is, Tiniuai — J2k^k- The magnetic field is 
gradually decreased by adjusting Tq and the critical point 
is attained at t = —Tq with B = 1. Then the geometric 
phase for the k th mode is 



1 



cos k + 1 



^/(cosfc + 1)2 



At the 

^critical 



(12) 

critical point the total geometric phase is 
— J2k^k{t = —Tq)- Finally, at t = 0, when the 
magnetic field is gradually turned off, the situation is a 
bit different. The configuration of the final state will de- 
pend on the number of kinks generated in the system 
due to phase transition at or near t — —Tq and as such 
it will depend on the quench time Tq [l9|. The num- 
ber of kinks is the number of quasi-particles excited at 
B — {) and is given hy H — X^fc-Pfc where pk, the ex- 
citation probability (for the slow transition) is given by 
the Landau Zener formula Pk ~ exp (— 27rrgfc^). As 
different pairs of quasi-particles {k,—k) evolve indepen- 
dently, for large values of Tq, it is likely that only one pair 
of quasi-particles with momenta (fco, —ko) will be excited 




FIG. 1: Color online, variation of geometric phase with time 
for different quenching time. The anisotropy parameter {a = 
0.5). Inset shows the same study but for q = (XX model). 

|l>-fe) 

where fco(= jj') corresponds to the minimum value of the 
energy A^. Thus the condition for adiabatic transition in 
a finite chain is given by, Tq >> Hence, well in the 
adiabatic regime, the final state at i = is given by 

\lpftnal >= |1 >k„ |0 >-ko 

n (cos ^|0 >k |0 >_fc -*e(2'*) Sin ^|1 >k |1 >^k) 

k,k^±ko 

This state is similar to direct product of only iV — 1 spins 
oriented along {2(j), Ok) where the state of the spin corre- 
sponding to momentum feg does not contribute to the ge- 
ometric phase. The total geometric phase of this state is 
given by F/i„a/(i = 0) = Y.k,k^±ko '^(l ~ cos^^) Now we 
study the geometric phase associated with the quench in- 
duced quantum phase transition. We study the variation 
of geometric phase (F^) and its derivative with respect 
to the quench induced magnetic filed {B) i.e. (^^) with 
time. We find the non-analytic behavior of the derivative 
ai t = Tq. The analytical expression for the derivative is 



dVkjt) 
dB 



\k) 



{{cOs{k)+t/Tqy 



\k)) 



3/2 



(13) 



In Fig. 1 we study the Berry phase F^ for different 
quench time Tg, for the XX and XY spin chain system. 
For the XX (inset) model, we find the step function like 
behaviour at t = Tq where the sharp transition occurs. It 
can be seen that the system shows quantum criticality at 
that point. The sharp transition disapperas for the finite 
anisotropy (XY) model. We observe from our study that 
the variation of F^ become flat for slower quenching rate. 
The upper panel of Fig. 2 shows the variation of F^ with 
a and t to get the whole picture of variation of F^. It 
reveals from our study that sharp transition occurs for 
a = only. The lower panel of the Fig. 2, shows the to- 
tal variation of The non-analytical behaviour for 
a — at t — —Tq helps us to predict XX criticality under 
the linear quenching process. The analyis with different 
values of Tq shows that the appearance of XX criticality 
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a t 



FIG. 2: Color online, the left planel of this figure is for the 
geometric phase and right one is for the derivative of the ge- 
ometric phase to study the quantum criticality. We present 
the vartion of these quantities with anisotropy parameter and 
time. 

is independent of r^, i.e., independent of fast and slow 
quenching rate. 

Finaly, we discuss an important aspect of the XX crit- 
icality. It is shown that the XX region of criticality is 
characterized by the existence of non-contractible geo- 
metric phase. The system is in the gapless excitations 
in the XX region of criticality (the gapless excitation is 
defined as bk'' operator). Here cos(^^) = B{t){Ako), 
fcoel, ...M. For finite Af, we can write Imia^o cos{9k) = 
±1 and ground state \g{a, B{t), cj)) >= i8)fc<fc„|0 >k 
|0 >-k '8)fc>fco|l >k |1 >-k- The ground state of the 
system acquires the geometric phase due to the vari- 
ation of (p for fixed a and B{t). \g > is the tensor 
product of AI qubit states. We observe from the Fig. 
1 and also from the analytical expression for F^; that 
limQ_>.o = 0, 27r. In the thermodynamic limit there is 
a solution, cos9k„ = 0, which leads to the result Tk = 
for every a > and hence limQ_>o Ffe^ — tt. There- 
fore by using the relation of total geometric phase, we 
can write limc^oliniAf^-oo^f- 7^ which shows the se- 
quence Ffc(a„)j^^jY is non-contractible in the thermody- 



namic limit. Hence, the non-contractible nature of the 
geometric phase associated with the XX criticality is 
proved. 

To conclude, we have studied the various quantum 
phases of the XY spin model and also the effect of 
linear quench on these quantum phases. The dynamics 
of the associated geometric phases are studied for 
slow and rapid quenching and an intimate connection 
between geometric phase and quantum criticality is 
established. We have predicted XX criticality and also 
non-contractible geometric phase at criticality. 
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